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We study the extremely neutron-rich nuclei for Z = 17 — 23, 37 — 40 and 60 — 64 regions 
of the periodic table by using axially deformed relativistic mean field formalism with 
NL3* paramctrization. Based on the analysis of binding energy, two neutron separation 
energy, quadrupole deformation and root mean square radii, we emphasized the speciality 
of these considered regions which are recently predicted islands of inversion. 



1. Introduction 

The main aim of theoretical models is to explain the available experimental results 
and predict the properties of the atomic nuclei through out the periodic table. 
A good description of the properties of known nuclei gives us more confidence in 
extrapolating to the yet unexplored areas of the nuclear chart. The two-neutron 
separation energy S^n systematically derived from the ground state binding energy 
(BE), reveal a new feature for the existence of islands of inversion in the e xotic 
neutron-rich regions of nuclear landscape. The Shell Model (SM) calculation DEE! 
is successful in nuclear structure theory. Although the application of this model in 
various regions explain the data quite well, it fails to reproduce the binding energy 
for some of the neutron-rich Ne, Na and Mg nuclei H Almost two decades ago 
Patra et al. SI performed the relativistic mean field (RMF) calculation with NLl 
parameter set and could explain the reason of failure of shell model for these nuclei. 
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One of their explanation is the large deformation of these nuclei which are not taken 
in the SM calculation. Rec ently apart from supporting the presently known islands 
around 31 Na ^ and 62 Ti ^ regions, the INM Model ^ predict one more region 
around Z = 60 of stability. It was suggested that these nuclei of Z = 17 - 23, N = 38 
- 42, Z= 37 - 40, N = 70 - 74 and Z = 60 - 64, N = 110 - 116 regions are deformed 
and form islands of inversion with more binding energy than their neighbouring 
family of isotopes. This prediction motivate us to study the properties of such nuclei 
and to investigate the possible reasons of the extra stability. In the present paper, 
we have done the calculation for these three regions by using the axially deformed 
RMF model. 



2. Theoretical Framework 

The relativistic mean field (RMF) model | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | become famous in 
recent years and have been applie d to finite nuclei and infinite nu clea r matter. We 
have taken the RMF Lagrangian 03 with NL3* parameter set in our study. 
This force parameter is successful in both /3-stable and drip-line nuclei. The La- 
grangian contained the term of interaction between meson and nucleon and also 
self-interaction of isoscalar-scalar sigma meson. The other mesons are isoscalar- 
vector omega and isovector vector rho mesons. The photon field is included to 
take care of coulombic interaction of protons. A definite set of coupled equations are 
obtained from the above Lagrangian are solved self-consistently in an axially de- 
formed harmonic oscillator basis. We start with the relativistic Lagrangian density 
for a nucleon-meson many-body system, 

C = WiV^ - M}ipi + ^&>ad»a - \m 2 y 

-■^.92ct 3 - ^53C 4 - g s ipiipiO- - ^tt^tt^ 

~B" V .B„ V + ^m 2 p &.R^ - gp^T^.R^ 

—F^F^ - e^ ^'^ ^V (1) 

Here sigma meson field is denoted by a, omega meson field by and rho meson 
field is denoted by p p . denotes the electromagnetic field, which couples to the 
protons, are the Dirac spinors for the nucleons, whose third components of isospin 
is T3 and g s , <?2, 33, 3w,C3, g p are the coupling constants. The center of mass (cm.) 
motion energy correction is estimated by the harmonic oscillator formula E c . m . = 
|(41^4 -1 / 3 ), where A is the mass number of the nucleus. The total quadrupole 
deformation parameter /?2 of the nucleus, can be obtained from the relation Q — 
Qn + Q P = \J ' ^-(^AR 2 ^), where Q n and Q p are the quadrupole moment for 
neutron and proton respectively and R is the nuclear radius. The root mean square 
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charge radius (Tv^), proton radius (r p ), neutron radius (r„) and matter radius (r m ) 
fl9l 

are given as lLMi . 



< rl >= i J p n (r±,z)rldT n , 



r c h 



< r >= 



r 2 v + 0.64, 



p(r^,z)r 2 dr, 



(2) 
(3) 
(4) 
(5) 



here all terms have own usual meaning. The total binding energy and ot her observ- 
ables are also obtained by using the standard relations, given in I 17 | 19 [ 



2.1. Pairing Correlation 

Pairing correlation is playing very crucial role in open shell nuclei. In our calculation 
we are using the Bardeen-Cooper-Schrieffer (BCS) pairing for determining the bulk 
properties like binding energy (BE), quadrupole deformation parameter and nuclear 
radii. The pairing energy can be given as: 

-i 2 



E,i 



-G 



i>0 



(6) 



where G is p airing for ce constant and vf , iq = 1 — vf are occupation probabilities 
190191 1991 

respectively 1 1 1 The simple form of BCS equation can be derived from the 
variational method with respect to the occupation number vf: 

ItiUiVi - A(uf - vf) = 0, (7) 

using A = G^2 i>Q UiVi. The above equation [7] is known as BCS equation for 
pairing energy. 

The occupation number is defined as: 



1 



n, = v. 



1 



A 



(8) 



- A) 2 + 

In our calculation we are dealing with the nuclei far away from beta stability 
line, so the constant gap for proton and neutron used here is valid in considered 
region . These constant gap equation for proton and neutron is taken from Ref. 
| 23 | 24 | wn j c h j s given as: 

A p = RB s e sI ~ tl2 /Z 1 ^ (9) 

and 

An^RB^-^/A 1 / 3 , (10) 
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with i?=5.72, s=0.118, t= 8.12, B s =l, and I = (N - Z)/(N + Z). 

In our present calculation, we have taken the constant pairing gap for all states 

| a >=| nljm > near the Fermi surface for the shake ol simplicity. As we know, 

if we go ne ar the very neutron drip line, then coupling to the continuum become 

important | 25 | 26| j n cage we snou u use the Relativistic Hartree-Bogoliubov 

(RHB) approach which is more accurate formalism for this region. But by using 

BCS pairing correlation model, it has been shown that the results from relativistic 

mean field BC S (RMF-BCS) approach is almost similar with the RHB formalism 
I27I28I29I30I3T1 

2.2. Choosing the Basis 

We divide our calculation between three regions: first region having Z—17 -23, and 
second region Z= 37-40 and third region Z—60 - 64. For checking the proper basis, 
we calculate the physical observables like binding energy, root mean square (rms) 
radii and quadrupole deformation parameter^)- So we have taken heavier nuclei 
from each region for example e9 V from region I, 119 Zr from region II and 1S5 Gd 
from region III. We have presented our calculations in table 12. 2\ with Nf^Nb 
=6 to 16, in the interval of 2, at the initial deformation of fa=0.2, using the NL3* 
parameter set. For 69 V nuclei; BE, rms radii and fa are almost same for Np=Ns 
> 10. It means, we can take N max — 10 for boson and Fermion harmonic basis for 
region I. In 119 Zr nucleus; these physical observables change from Nf=Nb— 10 to 
12. but become constant after Nf=Nb= 12. Then if we combine these two region 
and take N max =12, then we have enough space for both region. One can easily see 
that for 185 Gd nucleus N max =12 is not sufficient model space. It needs more space 
for calculation. Therefore, in this paper, we have taken N max —12 for region I, II 

IQOI 

and N max =14 for region III. Similar calculation are found in Ref. 

The physical observables like binding energy, root mean square(rms) radii and 
quadrupole deformation parameter (fa) does not much change when blocking is 
applied I^ZI. We have given our calculated results in table 12.21 The RMF values of 
the observables without blocking remains same with the values of blocking. If we 
see the table 12. 2\ the difference between blocking and without blocking is nearly 
less than 1 MeV. 



3. Choosing Reference Frame 

While comparing our binding energy results with macro-microscopic (MM) ap- 
proach, some important points needed to be stated. It is a known feature in MM 
models that the order of accuracy varies from region to region ^ i n the N-Z plane. 
The degree of disagreement is unacceptably large even slightly away from the known 
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Table 1. Calculated binding energy BE(MeV), root mean square r c h, r n , r p , r m and quadrupolc 
deformation parameter 02 at different bosonic and Fcrmionic basis harmonic quanta. Root mean 
square radius are in fm. During this calculation we have taken initial deformation parameter 0q 
equal to 0.2. 



Nucleus 


Basis 


BE 


Teh 


r n 


r p 








6 


508.810 


3.733 


4.103 


3.646 


3.956 


0.160 




8 


525.829 


3.822 


4.236 


3.738 


4.077 


0.220 




10 


530.960 


3.827 


4.301 


3.742 


4.123 


0.237 




12 


531.547 


3.831 


4.330 


3.747 


4.145 


0.247 




14 


531.651 


3.832 


4.343 


3.747 


4.154 


0.249 




16 


531.779 


3.831 


4.355 


3.747 


4.162 


0.251 


119 Zr 


6 


816.970 


4.302 


4.810 


4.227 


4.622 


0.139 




8 


904.886 


4.500 


4.903 


4.428 


4.749 


0.084 




10 


919.504 


4.518 


4.986 


4.447 


4.811 


0.055 




12 


922.347 


4.525 


5.018 


4.454 


4.836 


0.011 




14 


922.768 


4.525 


5.025 


4.454 


4.841 


0.006 




16 


922.868 


4.524 


5.031 


4.453 


4.844 


0.005 


185 Gd 


6 


1057.418 


4.783 


5.507 


4.716 


5.247 


0.139 




8 


1341.916 


5.155 


5.520 


5.092 


5.376 


0.106 




10 


1395.818 


5.280 


5.661 


5.219 


5.512 


0.118 




12 


1406.994 


5.294 


5.721 


5.233 


5.557 


0.122 




14 


1408.427 


5.289 


5.737 


5.228 


5.566 


0.114 




16 


1408.222 


5.289 


5.741 


5.228 


5.569 


0.114 



domain (See Figs. 7-9, Ref. '"'and Fig. 1, Ref. *^). On the other hand a microscopic 
formalism based on nuclear Lagrangian/Hamiltonian predicts physical observables 
through out the known/unknown territory of the periodic chart equally well. The 
parameters of these models have been determined by fitting the experimental data 
of few well known nuclei only. It is to be noted that the prediction of nuclei even in 
the known region are treated in an equal footing with the unknown region. There- 
fore, similar predictive power can be expected in the actual unknown region. 

4. Calculations and Results 

Relativistic mean field model have given very good result in f3 stable nuclei of the 
nuclear landscape. In this work we are analyzing the exot ic neutron drip line nuclei 
by using RMF model with recent well known NL3* parameter set. We obtain 
matter radii, quadrupole deformation parameter and ground state binding energies 
of these exotic nuclei of Z = 17-23, 37-40 and Z = 60-64 regions. The calculated 
results, like binding energy, radii, quadrupole deformation are given in tables (3- 
5) and the results are discussed in figures 1-10. In upcoming subsections we have 
described these results in detail. 
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Table 2. Calculated ground state binding energy BE(MeV), root mean square f^fr.; 'f'ni 

and 

quadrupole deformation parameter 02 without blocking and with blocking. Root mean square 
radius are in fm. 



Nucleus 


Without Blocking 


With Blocking 




BE 


Teh 


r n 


Tn 




#2 


BE 


T c h 


Tn 


Tn 
I J 






51 Cl 


385.05 


3.49 


4 


3.39 


3.81 


-0.25 


384.46 


3.49 


4 


3.39 


3.81 


-0.25 


63 Cl 


390.19 


3.63 


4.51 


3.54 


4.27 


0.31 


389.97 


3.63 


4.51 


3.54 


4.27 


0.31 


5i Ar 


402.94 


3.5 


3.94 


3.41 


3.76 


-0.23 


402.48 


3.5 


3.94 


3.4 


3.76 


-0.23 


63 Ar 


416.23 


3.65 


4.41 


3.56 


4.19 


0.23 


416 


3.65 


4.4 


3.56 


4.18 


0.22 


53 K 


423.78 


3.5 


3.96 


3.41 


3.77 





423 


3.51 


3.96 


3.41 


3.77 


-0.01 




441.36 


3.65 


4.3 


3.56 


4.09 


0.12 


440.79 


3.65 


4.3 


3.56 


4.09 


0.12 


55 Co 


446.46 


3.56 


3.98 


3.47 


3.8 





445.91 


3.56 


3.98 


3.47 


3.8 


-0.03 


65 Ca 


465.8 


3.7 


4.34 


3.61 


4.13 


0.14 


465.55 


3.7 


4.34 


3.61 


4.13 


0.14 


56 Sc 


461.32 


3.6 


3.97 


3.51 


3.8 





460.78 


3.6 


3.97 


3.51 


3.8 





66 Sc 


487.47 


3.74 


4.33 


3.66 


4.12 


0.18 


486.79 


3.74 


4.32 


3.66 


4.12 


0.18 


57 Ti 


475.64 


3.64 


3.96 


3.55 


3.81 


-0.1 


475.08 


3.64 


3.96 


3.55 


3.81 


-0.1 


57y 


484.57 


3.68 


3.91 


3.59 


3.79 


0.17 


483.79 


3.68 


3.91 


3.59 


3.79 


0.17 


68y 


529.77 


3.82 


4.3 


3.74 


4.12 


0.23 


528.93 


3.82 


4.29 


3.74 


4.11 


0.23 


103 Rb 


832.98 


4.42 


4.81 


4.35 


4.65 


-0.27 


832.42 


4.41 


4.81 


4.34 


4.65 


-0.26 


lw Rb 


853.43 


4.43 


4.94 


4.35 


4.75 


-0.06 


852.67 


4.43 


4.94 


4.35 


4.75 


-0.07 


105y 


863.76 


4.44 


4.79 


4.37 


4.64 


-0.23 


863.17 


4.45 


4.79 


4.38 


4.64 


-0.24 


107 Zr 


882.85 


4.47 


4.81 


4.4 


4.66 


-0.23 


882.4 


4.48 


4.81 


4.4 


4.66 


-0.24 



4.1. Binding Energy 

Binding energy (BE) is precisely observed from experiment which is responsible for 
stability and structure of the nuclei. The maximum binding energy corresponds to 
the ground state for a given nucleus and a ll other solutions are intrinsic excited 
states. We have given RMF, FRDM IMISSI INM grouri d s t a te binding energy in 
second, third and fourth column of the tables (tables 3-5) respectively. In fifth 
column AEi which indicates the binding energy difference between FRDM and 
RMF i.e. BE(FRDM) - BE (RMF) and in sixth column AE 2 indicates the binding 
energy difference between INM and RMF i.e. BE(INM) - BE(RMF). Last two 
columns for quadrupole deformation parameter (QDP)(/32) of RMF and FRDM 
model respectively. In this subsection we are comparing our RMF binding energy 
(BE) with INM (BE) ^and well established FRDM (BE) E5 results. 

In fig. 1(a), in Z = 17 — 23 region, we have plotted the binding energy difference 
AE for CI isotopes. We get AE 2 is zero means RMF and INM binding energies are 
nearly same at lower mass region, but if we go further, difference will increase in 
middle part and at A = 58, 59 again it goes to nearly zero, but diverges at higher 
mass region. If we compare our result with FRDM, then we got AEi nearly zero at 
lower mass region, but it diverges at higher mass region. In Fig 1(b), in case of Ar 
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isotopes. The RMF BE is not consistent with INM at lower mass region but we get 
AE 2 nearly zero at middle region at A — 54 - 60, which again diverges at higher 
mass region. If we compare our results with FRDM, then we get AEi nearly zero at 
lower mass region at A — 51 - 56 then the difference increases at higher mass region 
at A = 56 - 61. In Fig 1(c), RMF binding energy is very close to INM in lower 
mass region at A = 52, 54 then AE 2 increases within A = 55 — 60 in the middle 
region. Further it is very close to INM binding energy. If we compare our results 
with FRDM, then AEi tends to zero at A= 53 then AE\ increases further. In Fig 
1(d), in case of Ca, RMF binding energy is very close to INM and FRDM binding 
energy at lower mass region at A = 53, 56. RMF binding energy diverges from 
both model (INM and FRDM) in the middle part A = 56 - 62 and then matches 
at higher mass region. In Fig 1(e), in case of Sc, we got AE\ and AE2 are zero in 
lower mass, whereas both diverge in the middle part A = 53 - 63 then it further 
moves to zero. In the Fig 1(f), AE\ and AE2 are following the same trend as Fig 
1(e) but it diverges in A = 58 - 62, and at higher region RMF BE matches with 
FRDM and INM predictions. 

The binding energy difference for Rb isotopes is given in Fig. 2(a), the AEi has 
a large value at lower mass A = 103 - 107, then it tends to zero in higher region but 
if we compare RMF with INM results, AE 2 increases in lower mass region and go 
to zero in middle region then diverges at higher mass A = 107 - 114 region. In Fig. 
2(b), we plotted the AEi and AE 2 for Sr isotopes. We got same trend but RMF 
results diverges from INM at higher mass region while it closes to FRDM. In lower 
mass region RMF results are not matching with INM and FRDM results. Energy 
difference AE for Y nuclei isotopes are given in Fig. 2(c), again RMF results are 
not consistent with INM and FRDM results at lower mass A = 105 - 108, but at 
higher mass region it matches with INM and FRDM results. Fig. 2(d), represent 
AEi and AE 2 for Zr isotopes, from figure it is clear that our RMF results are not 
matching with INM and FRDM results. 

In Fig. 3, we have given AE (binding energy difference) for region Z — 60 - 64 
nuclei. For Nd isotopes, RMF binding energy is not consistent with FRDM at A = 
166 -180. Later on RMF binding energy is close to FRDM result for few isotopes 
A = 179 - 181 then again diverges at A = 182. When we compare our result with 
INM binding energy, the binding energy of RMF is very close to it at A = 168 - 170 
and later on diverges with increase in mass number. In case of Pm nuclei isotopes, 
which is plotted in Fig. 3(b), RMF result is not consistent with FRDM for the 
whole region. If we compare our result with INM region, it is consistent till A = 
168 - 172 and then diverges. In Sm isotopes, RMF binding energy is not consistent 
with FRDM in whole region. AE for Sm has followed the same trend as AE of 
Pm nuclei isotopes i.e. matches at lower mass region and diverge at higher mass 
region for both FRDM and INM results. In Eu nuclei isotopes, RMF BE does not 
consistent with FRDM and INM binding energy, in both case AE increases with 
mass numbers. 
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Fig. 1. Difference between the binding energies using RMF, Finite Range Droplet Model (FRDM), 
Infinite Nuclear Matter (INM) model (a) The circle represent for AEi (FRDM - RMF) (b) The 
square represent for AE2 (INM - RMF) model for different mass values of Z = 17 - 23 region. 
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Fig. 2. Same as Fig. 1 for Z = 37 - 40 region. 
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Fig. 3. Same as Fig. 1 for Z = 60 - 64 region. 



4.2. Quadrupole Deformation 

Quadrupole deformation parameter (QDP) 02 is directly connected to the shape 
of the nucleus. It is very common to say that if we go towards drip-line nuc lei, 
deformation will gradually increase but recently experimental paper of Tshoo 
explain that 22 is prolate in shape but 24 is spherical in structure. Keeping this 
result in our mind we have calculated the QDP 02 for recently predicted island of 
inversion region in nuclear landscape. Because of the unavailability of experimental 
data of these nuclei, we have compared our calculated QDP 02 with well stabilized 
FRDM ESkata. In Fig. 4, we have plotted the quadrupole deformation parameter 
fa for RMF and FRDM models as a function of mass number for Z = 17 - 23 
region. In CI case, QDP 02 continuously increases with the mass number, as shown 
in Fig. 4(a). In lower mass region CI isotopes are oblate and in higher mass region 
these are prolate and middle case A = 56 - 58, there is continuous shape change 
from oblate to prolate. If we compare RMF results with FRDM predictions then 
we get totally different result in FRDM. In FRDM, shape is suddenly changed from 
oblate to prolate (A = 54 ) and prolate to oblate (A = 57). Most of the CI isotopes 
are oblate in FRDM model. There are continuous changes in deformation but there 
is very small amount of energy difference (1 MeV) between ground state and first 
excited state. So we can say that other shape is also possible, But here we are taking 
only the ground state and neglecting the other possibility of shapes. In Ar case, 
most of the isotopes are oblate in lower mass region A = 52 - 57, and some are 
spherical at A = 59 - 60 then in higher region it again changes its shape from oblate 
to prolate in RMF model. When we compare with FRDM data, RMF is very close 
to FRDM except middle and high region in Fig. 4(b). FRDM is completely oblate 
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in shape over the region. In Fig. 4(c), we have plotted the QDP f3 2 for K isotopes. 
From figure it is very clear that most of the isotopes are spherical in shape. When 
we compare with FRDM data, it shows the same trend as RMF at A = 54-57 i.e. 
spherical in shape. In Ca, Sc, Ti case all are spherical in shape over all isotopes. 
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Fig. 4. Quadrupolc Deformation Parameter obtained from RMF(NL3*) (circle) compared with 
the FRDM(squarc) results for different isotopes of Z = 17 - 23 region. 



Deformation parameter for Rb isotopes are given in fig. 5(a). From the figure it 
is clear that most of the isotopes are spherical shape but in lower mass region A = 
103 - 105, it is oblate. If we see the QDP /?2 for Rb isotopes in FRDM model, we 
found that most are in prolate shape in lower mass region A = 103 - 110 then shape 
changes to oblate which is totally different from RMF result. Sr isotopes are given 
in fig. 5(b), for Sr isotopes, most of the nuclei are in spherical but in lower mass A 
= 104 - 106 are prolate and at A = 103 shape changes from prolate to oblate. If we 
see the result of FRDM, most of the Sr isotopes are prolate and in higher region it 
is spherical. RMF matches to FRDM at A = 104 - 106 and in higher region. Again 
we are getting spherical shape for A = 108 - 117 for Y isotopes in fig. 5(c). RMF 
matches only at A = 106, 107 and in higher mass A = 114 - 116. In Zr isotopes, 
It is spherical in shape at A = 109- 120 except A = 114, 115 as shown in the fig. 
5(d). If we go from A = 107 to 109, then we got a sharp shape change at A = 108 
i.e. oblate to prolate and again prolate to spherical. FRDM have prolate shape in 
lower mass region A = 107 - 113 and then changes to oblate in A = 114 - 120. 

In Fig. 6(a), for Nd isotopes, at A = 167 -174, both RMF and FRDM are prolate 
in shape but when we go further RMF change its shape to oblate in A = 175 - 179 
region while FRDM remains prolate in shape. In higher mass region A = 180 - 
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Fig. 5. Same as Fig. 4 for Z = 37 - 40 region. 



182 RMF goes oblate to nearly spherical and FRDM goes from prolate to oblate 
it means these two model are not consistent in A = 175 - 182. Isotopes of Pm are 
given in fig. 6(b). Here we get consistent result in RMF and FRDM model at A = 
167 - 175. Then RMF changes to oblate in higher mass region and again change to 
nearly spherical while FRDM does not change the shape. In Sm case, both models 
match to each other in A = 168 - 176 then RMF goes to oblate and spherical shape 
where FRDM does not follow the RMF trend except A = 181, 182. In Eu isotopes 
as shown in fig.6(d), both RMF and FRDM show consistency at A = 169-178, later 
on RMF changes to oblate at A= 179-182. FRDM is matching with RMF at A 
=181,182 only at higher mass region. 



4.3. Nuclear Radius 

In this subsection we are concentrating on the neutron radius (r n ), proton radius 
(r p ), charge radius (r c h) and matter radius (r m ) which are calculated by using 
RMF(NL3*) formalism. In Fig. 7, we have plotted the r n , r p , r c h, and r m for 
CI, Ar, K, Ca, Sc, Ti and V nuclei. In Z = 17 - 23 region, all the radii increase 
monotonically with mass number. In Fig. 8, we have plotted the r n , r p , r c h, r m 
with mass number for Z = 37-40 region. In Rb isotopes, there is a sharp fall in 
radii till A=106, then radii increase monotonically. In Sr isotopes, the radii follow 
same trend as Rb isotopes but in this case fall at A=107, then the radii increase 
monotonically. In Y isotopes, the radii follow a jump at A=106 and remain constant 
upto A = 107, then decrease at A=108. Later on the radii follow the same trend 
means the radii increase monotonically with mass number. In Zr isotopes, the radii 
increase and it follow a jerk at A = 108 then go down at A = 109 and later on 
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A 

Fig. 6. Same as Fig. 4 for Z = 60 - 64 region. 

increase. In Fig. 9, we have plotted the radius curve for Z = 60 - 64 region, in the 
case of Nd, Pm isotopes the radii increase monotonically with atomic number. In 
Sm (Fig. 9c) isotopes, we get a small jerk in A = 176 while in Eu (Fig. 9d) isotopes 
this jerk arises at A = 177-178, then increases monotonically saying a change in 
the deformation of the nuclei. 

4.4. Two-neutron separation energy 

The two-neutron separation energy S^^NjZ) = BE(N, Z) - BE(N-2, Z) is shown in 
Fig. 10. The S2n values decrease gradually with increase in neutron number. It is 
indeed satisfying to note that in th e recent years strong evidences both experimental 
and theoretical have emerged supporting the existence of this island of inversion 
centering around 62 Ti. We can predict the stability of these nuclei by S^n energy. If 
S2n is large, it means nuclei will be stable with two-neutron separation. As shown 
in first part (a) of Fig. 10, in Z = 17-23 region, we are getting a sharp down curve 
for all the members of this region at N = 42. So we can say that this may be the 
neutron magic number in this neutron-dripline nuclei. In S^n plot for Ti, we are 
getting a small considerable jerk at N = 44. This shows the extra stability of nuclei. 
In Sc, S2n plot follow the same trend as in Ti, but the magnitude is very small. In 
other cases, i.e. Z = 17, 18, 19, 20 region, there is no any local extra stability. In 
second part (b), we are getting a sharp down curve at N = 68 for all the members 
of this region. In Sr, Rb, there is a small jerk at N = 74. In other cases, there is no 
local stability. In third part (c), for Z = 60 - 64 region, there is a sharp fall at N = 
112 for all the members of this region. We get local extra stability in Nd, Pm and 
other nuclei also follow nearly the same trend. 
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Table 3. We have given ground state binding energy (BE) for RMF(NL3*), FRDM, INM and 
comparison between calculated BE using RMF(NL3*)model with FRDM(AEi) and INM (A£ 2 ) 
and /?2 for Z = 17 - 23 region. Binding energies are in MeV. 



nucleus 


BE(RMF) 


BE(FRDM) 


BE(INM) 


AEi 


AE 2 


£ 2 (RMF) 


#2 (FRDM) 


&i C7 


385.1 


383.1 


387.2 


-2.02 


2.18 


-0.256 


-0.307 


52 CI 


386.9 


382.9 


386.5 


-4.07 


-0.44 


-0.221 


-0.357 


53 C7 


388.8 


384.3 


387.7 


-4.57 


-1.13 


-0.177 


-0.389 


5i Cl 


390.4 


383.8 


385.9 


-6.56 


-4.48 


-0.155 


-0.406 


55 CI 


391.9 


384.7 


388.7 


-7.15 


-3.20 


-0.134 


0.452 


56 CI 


393.0 


384.0 


388.5 


-9.01 


-4.59 


-0.085 


0.454 


57 Cl 


394.7 


385.1 


391.9 


-9.64 


-2.80 


0.001 


0.456 


58 CI 


393.5 


385.4 


391.8 


-8.04 


-1.70 


0.071 


-0.262 


59 CI 


392.7 


386.6 


394.6 


-6.17 


1.83 


0.160 


-0.302 


60 CI 


392.2 


385.2 


394.3 


-6.99 


2.12 


0.216 


-0.301 


ei Cl 


391.8 


386.6 


396.8 


-5.14 


5.03 


0.254 


-0.294 


62 CI 


390.9 


384.5 


395.8 


-6.44 


4.86 


0.283 


-0.294 


63 CI 


390.2 


383.5 


397.4 


-6.67 


7.17 


0.310 


-0.297 


49 Ar 


395.9 


396.3 


399.9 


0.44 


4.01 


-0.201 


-0.223 


50 Aj , 


399.5 


401.0 


405.3 


1.49 


5.84 


-0.217 


-0.248 


51 Ar 


399.5 


402.0 


406.1 


2.57 


6.59 


-0.214 


-0.281 


52 Ar 


402.9 


405.9 


408.7 


2.98 


5.74 


-0.235 


-0.306 


53 Ar 


406.2 


406.1 


408.5 


-0.09 


2.22 


-0.247 


-0.339 


54 Ar 


408.7 


408.9 


410.4 


0.2 


1.72 


-0.219 


-0.357 


55 Ar 


411.1 


409.6 


409.0 


-1.47 


-2.05 


-0.182 


-0.307 


56 Ar 


413.3 


412.9 


411.8 


-0.36 


-1.46 


-0.159 


-0.237 


57 Ar 


415.4 


411.9 


412.2 


-3.5 


-3.18 


-0.137 


-0.229 


58 Ar 


416.9 


413.6 


415.6 


-3.36 


-1.37 


-0.084 


-0.255 


59 Ar 


419.1 


412.8 


415.8 


-6.33 


-3.38 


0.000 


-0.271 


60 Ar 


418.2 


414.6 


419.3 


-3.55 


1.16 


0.001 


-0.285 


61 Ar 


417.9 


413.9 


419.2 


-3.98 


1.36 


-0.105 


-0.280 


62 Ar 


416.8 


416.0 


422.3 


-0.74 


5.56 


-0.132 


-0.292 


63 Ar 


416.3 


413.8 


422.1 


-2.42 


5.89 


0.211 


-0.294 


53 R 


423.8 


422.5 


425.0 


-1.33 


1.22 


0.001 


-0.323 


54 K 


427.1 


423.7 


425.1 


-3.41 


-1.98 


-0.079 


-0.373 


55 K 


430.4 


426.5 


427.8 


-3.89 


-2.64 


-0.102 


-0.398 


56 K 


433.5 


428.2 


427.8 


-5.30 


-5.72 


-0.094 


-0.138 


57 R 


436.4 


431.1 


431.5 


-5.32 


-4.93 


-0.063 


-0.129 


58 K 


439.5 


431.4 


431.8 


-8.19 


-7.73 


0.000 


-0.155 


59 R 


442.8 


434.3 


436.0 


-8.41 


-6.74 


0.000 


-0.008 


60 R 


442.4 


434.5 


436.6 


-7.91 


-5.80 


0.000 


-0.272 


61 K 


441.6 


436.7 


440.6 


-4.88 


-1.03 


-0.045 


-0.283 


62 K 


441.4 


436.6 


441.2 


-4.76 


-0.20 


0.133 


-0.298 




441.4 


438.0 


444.7 


-3.31 


3.37 


0.173 


-0.297 


54 Ca 


442.7 


442.7 


443.7 


-0.07 


1.07 


0.000 


0.000 


55 Ca 


446.5 


444.5 


445.3 


-1.93 


-1.19 


0.000 


0.000 


56 Ca 


450.2 


449.2 


448.8 


-1.00 


-1.34 


0.000 


-0.018 


57 Ca 


454.0 


449.7 


450.0 


-4.26 


-3.95 


0.000 


-0.070 


58 Ca 


457.9 


454.5 


453.1 


-3.40 


-4.79 


0.000 


-0.007 


59 Ca 


461.9 


455.0 


453.9 


-6.92 


-7.95 


0.000 


-0.007 


60 Ca 


465.6 


458.2 


457.9 


-7.44 


-7.68 


0.000 


0.000 


61 Ca 


465.9 


458.9 


459.0 


-7.02 


-6.93 


0.003 


-0.014 


62 Ca 


465.6 


461.7 


463.3 


-3.93 


-2.36 


0.000 


0.008 


63 Ca 


465.3 


460.9 


464.3 


-4.4 


-1.02 


0.019 


0.045 


64 Ca 


465.6 


463.4 


468.4 


-2.21 


2.76 


0.126 


0.045 


65 Ca 


465.8 


462.6 


469.2 


-3.21 


3.39 


0.154 


0.071 


66 Ca 


466.0 


464.6 


472.6 


-1.33 


6.64 


0.172 


0.062 



continued 
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Fig. 7. The charge radii r c h (circle), The neutron radius r n (square), The proton radius r p (dia- 
mond), the rms radii r m of matter distribution (triangle up) for different isotopes of Z = 17 - 23 
region using the RMF(NL3*) formalism. 




5. Discussion and Remarks 

Taking RMF as a reference, we evaluate AEi and Ai?2- Analyzing fig. 1, we find 
both AEi and AE2 similar for all the considered six nuclei, CI to Ti. The large value 



25, 2013 2:11 WSPC/INSTRUCTION FILE island 



Island of inversion.... 15 



5.8 

5.6 

5.4 

^5.2 

£ 5 
^5.8 

5.6 

5.4 

5.2 



I I I I | I I I I | I I I I | I I I I 

:(a) 




; ^«-» «>♦<»<♦■» 
I I I I | I I I I | I I I I | I I I I 



(c) 



(d) 



Sm 



_i i i i i i i i i i i 



t — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 

"(b) " 



Pm 



I I I I | I I I I | I I I I | I I I I 




I ■ ■ ■ ■ I ■ ■ ■ ■ I 



170 175 180 185 170 175 180 185 

A 



Fig. 9. Same as Fig. 7 for Z = 60 - 64 region 




Fig. 10. The two-neutron separation energy S2 n for different isotopes of Z = 17 - 23, 37 - 40, 60-64 
region by using RMF(NL3*) formalism. 



of AEi and AE 2 at middle of the region shows the speciality of these nuclei, except 
CI isotopes [fig. 1(a)]. All other isotopes show similar trend with INM and FRDM. 
From fig. 2 and fig. 3, AEi are almost constant, if one extends the calculation to 
higher mass number in isotopic chain. On the other hand, calculated AE 2 goes on 
increasing with A. In this situation, the predictive power of RMF, FRDM and INM 
are questionable. For example, (1) if we consider RMF as the absolute reference 
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continued 



nucleus 




T3T7 /'T?T>T^A R\ 

Drj(r KDJV1 ) 


T> T? j T1\T1\ l\ \ 

r>Hj(llNM j 






a /pi\ 'XT? \ 


P2(r HUM j 


55 c 

be 


A C 7 1 

457.1 


A C O 7 

456.7 


A C 7 C 

457.5 


-0.40 


A A O 

0.43 


A AAA 

0.000 


A A70 

-0.078 


56 c 

be 


4ol.o 


A C A A 

459.4 


A C A A 

459.0 


1 A 1 

-1.91 


o o o 

-2.28 


A A A 1 

0.001 


A 1 A A 

-0.104 


57 c 

DC 


a otz a 

465.6 


AO A 1 

464.1 


A OO 7 

463.7 


1 CO 

-1.53 


-1.96 


A A A 1 

0.001 


A C\C\0 

-0.096 


58 o„ 

be 


/i 7 a a 
4^0.0 


A OO O 

466.3 


/l C A 

465.0 


O O C 

-3.65 


A A1 

-4.91 


A AAA 

0.000 


A 1 AC 

-0.105 


59 c 

be 


474.4 


A 7 A O 

470.8 


i a c\ o 

469.6 


O O^t 

-3.62 


A O O 

-4.83 


A AAA 

0.000 


A A7A 

-0.079 


60 c 

be 


/i to n 

478.9 


A TO 7 

472.7 


A 71 A 

471. (J 


O O /I 

-6.24 


7 07 

-7.87 


A AAA 

0.000 


A A A O 

-0.042 


61 c„ 

be 


A O O 1 

483.1 


/I 77 A 

4f f\0 


/I 7C O 

4^5.8 


-6.13 


7 

-7.2b 


U.000 


A Al C 

-0.015 


62 c 

be 


484.1 


/I 77 O 

477.3 


/i 77 C 

477.5 


C 77 

-6.77 


O CO 

-6.58 


0.000 


A A A 7 

-0.047 


63 c 
DC 


A O A O 

484.3 


A OA /I 

480.4 


A OO O 

482.3 


O A /I 

-3.94 


O A C 

-2.05 


A AAA 

0.000 


A A C O 

0.053 


64 o„ 

be 


/IOC /I 

485.4 


A O 1 O 

481.2 


AOO O 

483.6 


A 1 O 

-4.18 


1 71 

-1.71 


A 1 AO 

0.108 


A A70 

0.0<3 


65 c 
DC 


/l oo c 

486.6 


/I O O O 

483.8 


A OO A 

488.0 


O O 1 

-2.81 


1 1 

1.31 


A 1 OO 

0.166 


A AAA 

0.090 


66 o„ 
DC 


/I 07 C 

487.5 


/I O A 1 

484.1 


A OA O 

489.2 


o o c 
-3.35 


1 7A 

-1.70 


A 1 O 1 

0.181 


A 1 1 O 

0.11b 




ioc c 

435.5 


/I O O 7 

438.7 


/l 07 

43 f .6 


O OA 

3.20 


O 1 A 

2.19 


A AAC 

0.005 


A AAA 
0.000 


51 TV 


A A 1 7 

442.7 


/I /I /I o 

444.8 


/l A A A 

444.0 


O A7 

2.0 7 


1 o o 

1.28 


A A C\0 

0.006 


A AAA 
0.000 


52ti • 
1 I 


/l /I O 7 

448.7 


A CO 7 

452.7 


/I c o 

452.6 


O A7 

3.97 


o o o 
3.83 


A AAA 

0.000 


A AAA 
0.000 


53 ta - 
i I 


/l C /I c 

454.5 


/I C 7 1 

457.1 


/I C 7 A 

457.4 


O OO 

2.66 


o o o 
2.88 


A AA1 

0.001 


A AAA 
0.000 


54t^ ■ 
I J 


A O A A 

460.0 


AO A O 

464.2 


A O A C 

464.5 


A 1 O 

4.12 


A A C 

4.45 


A A A O 

0.048 


A AAA 

0.000 


55 ta 
1 I 


465.4 


/I £; 7 A 

46 f .9 


A OO 1 

468.1 


O /I 7 

2.47 


O OO 

2.b3 


A AO C 

0.085 


A 1 O A 

0.134 


56^ ■ 
1 I 


A 7A 7 

470.7 


/I 7 /I /I 

474.4 


/I 7 A 1 

474.1 


O A 

3.69 


O A A 

3.40 


A 1 AO 

0.103 


A 1 O C 

0.135 


57ta 


4^5.7 


/I 77 1 

47 *M 


/I 7/2 7 

4^0.7 


1 /I c 

1.45 


1 c\o 

1.0b 


a 1 n a 
-0.104 


A 1 O C 

0.135 


58t^ • 
1 I 


/i OA 

480.6 


/I O O 1 

483.1 


/i OO A 

482.0 


O A 7 

2.47 


1 A A 

1.44 


A AAC 

-0.095 


A 1 A C 

-0.105 


59 ta 


/IOC o 

485.3 


/IOC c 

485.5 


A O A 1 

484.1 


A 1 OO 

0.122 


1 OA 

-1.29 


A A70 

-0.0 (O 


A 1 AC 

-0.105 


60t^ • 


/|AA (\ 

490.0 


/I A 1 A 

491. U 


A O O A 

488.9 


1 A7 

1.07 


A AO 

-0.98 


A AA1 
0.001 


A A7A 

-0.079 


61 TV 


/l A C A 

495. (J 


A AO 1 

493.1 


/l A1 A 

491.0 


1 07 

-1.87 


A AC 

-4.05 


A AAA 
0.000 


A A1 O 

-0.018 




/l AA 7 

499.7 


A A7 A 

497.9 


495.9 


1 7C 

-1.75 


O 77 

-3.77 


A AAA 

0.000 


A AAA 

0.000 




a 




AQR 1 




-o.ou 


ft ftftft 

u.uuu 


ft ft AO 


64 Ti 


502.4 


503.4 


503.2 


1.05 


0.86 


-0.035 


0.027 


65 Ti 


505.5 


503.8 


504.7 


-1.65 


-0.80 


0.167 


0.062 


57y 


484.6 


486.5 


486.7 


1.94 


2.09 


0.167 


0.181 


58 y 


490.0 


490.5 


490.1 


0.45 


0.13 


0.150 


0.163 


59 y 


495.2 


496.7 


496.1 


1.53 


0.95 


-0.125 


0.162 


60 y 


500.4 


499.8 


498.9 


-0.60 


-1.54 


-0.107 


-0.130 


61y 


505.4 


505.5 


504.4 


0.05 


-1.07 


-0.071 


-0.104 


62y 


510.8 


508.6 


506.8 


-2.17 


-3.95 


0.000 


-0.044 


63y 


515.9 


513.8 


512.1 


-2.11 


-3.80 


0.000 


0.018 


64 y 


518.3 


515.5 


514.6 


-2.89 


-3.79 


0.000 


0.053 


65y 


520.9 


520.1 


519.7 


-0.80 


-1.20 


0.108 


0.053 


66y 


524.5 


522.5 


521.8 


-1.98 


-2.76 


0.184 


0.155 


67y 


527.7 


525.9 


526.6 


-1.79 


-1.02 


0.220 


0.163 


68y 


529.8 


527.5 


528.3 


-2.33 


-1.53 


0.233 


0.161 


69y 


531.5 


530.9 


532.5 


-0.70 


0.98 


0.249 


0.169 



frame, then the large discrepancy of AE2 with mass number indicate the failure of 
INM near the drip- line region, or vice versa, (2) similarly if we analyze for AEi, 
it is somewhat constant with RMF for the entire region considered in the present 
paper. As we have discussed in the subsection 3, the RMF is based on microscopic 
origin in mesons and nucleons level. Except few fitted nuclei, all others masses, 
radii and quadrupole deformation are the predicted results for a large region of the 
periodic chart. The RMF results are found to be good for almost all the known 
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Table 4. Same as table 3 for Z = 37 - 40 region. 



nucleus 


BE(RMF) 


BE(FRDM) 


BE(lNM) 




A£ 2 


fh (KMF) 


P2 ("RDM) 


wi Rb 


833.0 


837.5 


836.9 


4.50 


-0.58 


-0.270 


0.341 


104 Rb 


836.3 


840.0 


838.9 


3.69 


2.53 


-0.276 


0.333 


105 Rb 


839.4 


844.2 


842.8 


4.78 


3.36 


-0.284 


0.334 


106 Rb 


843.0 


846.4 


844.5 


3.33 


1.51 


0.000 


0.335 


w7 Rb 


846.9 


849.7 


847.5 


2.81 


0.67 


0.000 


0.336 


10S Rb 


849.2 


851.4 


848.2 


2.24 


-1.00 


0.000 


0.356 


109 Rb 


851.0 


854.0 


850.9 


3.04 


-0.07 


0.015 


0.373 


110 Rb 


853.3 


854.9 


851.4 


1.54 


-1.95 


-0.065 


0.287 


in Rb 


855.5 


856.4 


854.3 


0.92 


-1.17 


-0.08 


-0.333 


112 Rb 


857.5 


857.0 


854.6 


-0.47 


-2.87 


-0.083 


-0.340 


113 Rb 


859.6 


860.9 


857.2 


1.23 


-2.45 


0.097 


-0.145 


114 Rb 


861.7 


862.1 


856.9 


0.42 


-4.71 


0.103 


-0.135 


103 Sr 


844.9 


850.1 


849.4 


5.19 


4.49 


-0.167 


0.360 


104 Sr 


848.7 


855.1 


854.5 


6.47 


5.88 


0.416 


0.361 


W5 Sr 


852.5 


857.8 


857.1 


5.27 


4.60 


0.417 


0.353 


ioe Sr 


856.2 


862.6 


861.4 


6.41 


5.26 


0.420 


0.354 


i07 Sr 


860.2 


864.8 


863.6 


4.58 


3.44 


0.000 


0.345 


ios Sr 


864.2 


868.8 


866.9 


4.53 


2.71 


0.000 


0.346 


109 Sr 


867.0 


870.4 


868.2 


3.40 


1.26 


0.000 


0.375 


uo Sr 


869.1 


873.5 


871.7 


4.32 


2.54 


0.000 


0.393 


m 5r 


871.7 


874.6 


873.0 


2.89 


1.36 


-0.058 


0.393 


112 Sr 


874.2 


876.9 


876.9 


2.80 


2.72 


-0.073 


-0.300 


113 Sr 


876.6 


878.9 


878.0 


2.28 


1.45 


-0.080 


-0.153 


iu Sr 


878.8 


882.3 


881.0 


3.49 


2.18 


-0.083 


-0.153 


115 Sr 


880.9 


883.6 


882.1 


2.69 


1.15 


-0.082 


-0.144 


116 Sr 


882.9 


886.5 


883.5 


3.57 


0.51 


-0.078 


-0.136 


117 Sr 


884.9 


887.4 


883.0 


2.51 


-1.90 


0.075 


-0.128 


118 g r 


886.8 


890.3 


879.9 


3.45 


-6.88 


0.059 


-0.120 


119 Sr 


888.68 


892.6 


878.9 


3.89 


-9.75 


0.001 


-0.008 


120 Sr 


890.6 


895.2 


880.4 


4.54 


-10.28 


0.000 


0.000 


105y 


863.8 


868.9 


868.9 


5.10 


5.13 


-0.230 


0.372 


106y 


867.9 


872.2 


871.9 


4.32 


4.09 


0.419 


0.364 


107y 


871.9 


877.0 


876.3 


5.11 


4.37 


0.421 


0.364 


108y 


876.3 


879.9 


878.6 


3.66 


2.31 


0.000 


0.356 


109y 


880.7 


884.1 


882.6 


3.40 


1.88 


0.000 


0.357 


llOy 


883.8 


886.3 


884.3 


2.46 


0.53 


0.000 


0.386 


lily 


886.4 


889.6 


888.2 


3.24 


1.86 


0.000 


0.394 


112y 


888.8 


890.9 


890.1 


2.09 


1.23 


0.000 


0.404 


113y 


891.4 


894.3 


893.9 


2.91 


2.51 


-0.047 


0.384 


114y 


894.1 


896.1 


895.4 


1.97 


1.34 


-0.071 


-0.153 


115y 


897.9 


899.6 


898.7 


1.73 


0.84 


-0.084 


-0.153 


116y 


899.2 


901.7 


899.2 


2.47 


0.03 


-0.084 


-0.153 


117y 


901.5 


904.7 


900.6 


3.15 


-0.941 


0.064 


-0.145 


107 Zr 


882.9 


887.9 


887.5 


5.04 


4.60 


-0.231 


0.364 


108 Zr 


886.9 


893.4 


892.4 


6.50 


5.50 


0.417 


0.365 


109 Zr 


891.6 


896.5 


895.2 


4.87 


3.59 


0.000 


0.357 


no Zr 


896.3 


901.2 


899.8 


4.93 


3.59 


0.000 


0.368 


lll Zr 


899.8 


903.5 


901.6 


3.70 


1.86 


0.000 


0.387 


U2 Zr 


902.7 


907.6 


906.4 


4.83 


3.65 


0.000 


0.395 


113 Zr 


905.5 


908.9 


908.4 


3.37 


2.87 


0.000 


0.404 


m Zr 


908.5 


913.3 


912.8 


4.79 


4.35 


-0.167 


-0.197 


115 Zj , 


911.2 


915.4 


914.3 


4.19 


3.03 


-0.154 


-0.170 


116 Zr 


914.1 


919.5 


917.8 


5.39 


3.66 


-0.098 


-0.162 


117 Zr 


916.9 


921.6 


918.9 


4.72 


2.11 


-0.097 


-0.153 


118 Zr 


919.6 


925.3 


920.8 


5.67 


1.12 


0.047 


-0.153 


iw Zr 


922.5 


926.9 


920.7 


4.49 


-1.81 


0.001 


-0.146 
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Table 5. Same as table-3 for Z = 60 - 64 region. 



nucleus 


BE(kMF) 


BE(FKUM) 


BE(lNM) 


A Ei 


AA 2 


02 (RMF) 




w 'Nd 


1310.65 


1317.3 


1314.9 


6.65 


4.25 


0.31 


0.28 


16S Nd 


1313.91 


1321.1 


1315.3 


7.19 


1.39 


0.3 


0.28 


im Nd 


1317.15 


1323.2 


1315.5 


6.06 


-1.64 


0.29 


0.28 


170 Nd 


1320.27 


1326.9 


1318.4 


6.63 


-1.87 


0.29 


0.27 


171 Nd 


1323.12 


1328.9 


1319.5 


5.79 


-3.62 


0.28 


0.26 


172 Nd 


1325.76 


1332.5 


1322.5 


6.74 


-3.26 


0.28 


0.25 


173 Nd 


1328.16 


1334.4 


1323.4 


6.24 


-4.76 


0.26 


0.24 


174 Nd 


1330.42 


1337.8 


1325.9 


7.38 


-4.52 


0.25 


0.22 


17b Nd 


1333.03 


1339.4 


1326.6 


6.37 


-6.43 


-0.19 


0.21 


176 Nd 


1335.51 


1342.6 


1328.6 


7.09 


-6.91 


-0.18 


0.19 


177 Nd 


1338.14 


1343.9 


1328.5 


5.76 


-9.64 


-0.13 


0.17 


178 Nd 


1341.16 


1347 


1329.9 


5.85 


-11.25 


-0.12 


0.16 


179 Nd 


1344.15 


1347.6 


1329.5 


3.45 


-14.65 


-0.11 


0.18 


iso Nd 


1347.07 


1350.6 


1330.6 


3.53 


-16.47 


0.09 


0.17 


181 Nd 


1350.23 


1352 


1329.9 


1.77 


-20.33 


0.08 


0.15 


182 Nd 


1353.25 


1354.9 


1331.1 


1.65 


-22.15 


0.07 


0.11 


167 Pm 


1321.05 


1326.8 


1323.1 


5.75 


2.05 


0.33 


0.28 


168p m 


1324.48 


1329.6 


1324.3 


5.12 


-0.18 


0.32 


0.28 


169p m 


1327.84 


1333.5 


1327.9 


5.66 


0.06 


0.31 


0.28 


170p m 


1331.15 


1336.1 


1329.4 


4.95 


-1.75 


0.31 


0.28 


171 Pm 


1334.34 


1339.8 


1332.8 


5.46 


-1.54 


0.3 


0.28 


172 Pm 


1337.35 


1342.2 


1334.3 


4.86 


-3.05 


0.29 


0.27 


173 Pm 


1340.15 


1345.7 


1337.3 


5.56 


-2.85 


0.29 


0.26 


17ip m 


1342.6 


1348 


1338.5 


5.4 


-4.1 


0.28 


0.25 


175 Pm 


1345.01 


1351.5 


1341.5 


6.49 


-3.51 


-0.2 


0.22 


176p m 


1347.89 


1353.5 


1342.4 


5.61 


-5.49 


-0.2 


0.21 


177 Pm 


1350.6 


1356.7 


1344.8 


6.1 


-5.8 


-0.19 


0.2 


178 Pm 


1353.15 


1358.5 


1345.4 


5.35 


-7.75 


-0.18 


0.17 


179 Pm 


1355.95 


1361.6 


1347.3 


5.65 


-8.65 


-0.13 


0.16 


1S0 Pm 


1359 


1362.6 


1346.9 


3.61 


-12.09 


-0.12 


-0.19 


181p m 


1362.1 


1365.7 


1348.5 


3.6 


-13.6 


0.1 


-0.18 


le8 Sm 


1334.23 


1340.5 


1337.5 


6.27 


3.27 


0.34 


0.28 


169 Sm 


1337.85 


1343.3 


1339.6 


5.45 


1.75 


0.33 


0.28 


170 Sm 


1341.32 


1347.8 


1343.5 


6.48 


2.18 


0.32 


0.28 


171 Sm 


1344.71 


1350.4 


1345.5 


5.69 


0.79 


0.31 


0.28 


172 Sm 


1348.05 


1354.6 


1349.3 


6.55 


1.25 


0.31 


0.28 


173 Sm 


1351.24 


1357 


1350.8 


5.77 


-0.43 


0.3 


0.27 


17i Sm 


1354.19 


1361.1 


1354.1 


6.91 


-0.09 


0.3 


0.26 


175 Sm 


1356.77 


1363.4 


1355.5 


6.63 


-1.27 


0.29 


0.25 


176 Sm 


1359.59 


1367.3 


1358.5 


7.71 


-1.09 


0.35 


0.23 


177 Sm 


1362.3 


1369.4 


1359.7 


7.11 


-2.6 


-0.2 


0.22 


178 Sm 


1365.2 


1372.9 


1362.3 


7.7 


-2.9 


-0.2 


0.2 


179 Sm 


1368 


1374.6 


1362.9 


6.6 


-5.1 


-0.19 


0.17 


180 Sm 


1370.67 


1378.1 


1365.1 


7.43 


-5.57 


-0.17 


0.16 


181 Sm 


1373.66 


1379.2 


1365.1 


5.55 


-8.56 


-0.13 


-0.19 


182 Sm 


1376.79 


1382.7 


1367.3 


5.91 


-9.49 


-0.12 


-0.18 


183 Sm 


1380.03 


1384.5 


1368.6 


4.47 


-11.43 


0.1 


-0.16 


169 Eu 


1346.83 


1351.2 


1349.6 


4.37 


2.77 


0.34 


0.28 


170 Eu 


1350.58 


1354.6 


1351.8 


4.02 


1.22 


0.34 


0.28 


171 Eu 


1354.2 


1359.1 


1356.4 


4.9 


2.2 


0.33 


0.28 


l72 Eu 


1357.72 


1362.2 


1358.8 


4.48 


1.08 


0.32 


0.28 



continued 
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continued 



nucleus 


Drj(ixMr ) 


T3T7 /'T?T>T^A i\\ 
JDX^r KlJlVi ) 


T> T?' I T1\T1\ l\ \ 

r>rj(HNM ) 


A TP 


A TP 


O /T> 1\ TT? \ 
p 2 (HM1 J 


Q 1 T?T> TM\ T\ 
P2\r KJJM J 


173 TP 

hu 


1 of 1 oo 

1361.23 


1 Off A 


1 0/""*O A 

1362.9 


5.17 


1.67 


A O 1 

0.31 


A OO 

0.28 


174 tp 

hu 


1364. 65 


1 OfC\ o 

1369.2 


1365.1 


A 

4.56 


A A E 

0.45 


A O 1 

0.31 


A OO 

0.28 


175 TP 

hu 


1367.8 


1 O TO O 

1373.3 


1 Of O f 

1368.5 


5.5 


A O 

0.8 


A O 1 

(J. 31 


A OT 

0.27 


176 tp„. 
hu 


1 OTA £?0 

13*0.63 


1 OT/? 

137b 


1 OTA 1 

137U.1 


E OT 

5.37 


A C O 

-0.53 


A O 

0.3 


A 

0.2b 


177 

' hu 


1 070 71 

1373.71 


1 O TO A 

1379.9 


1 O TO O 

1373.3 


o 

6.2 


A A 

-0.4 


A O E 

0.35 


A O A 

0.24 


178 

hu 


1375.37 


1 O OO K. 

1382.5 


1 O T A T 

1374.7 


f 1 O 

6.13 


1 f T 

-1.67 


A Of 

0.36 


A OO 

0.22 


179 Z7„, 
till 


1 OTA A 1 

137y.41 


1 O O C A 

1385. y 


1 OTT T 

13 (7.7 


/? /i a 

0.49 


1 T1 

-1.71 


A O 

-0.2 


A 01 

0.21 


180 

hu 


1 o oo o o 

1382.38 


1 o o o o 

1388.2 


1 O TO A 

1378.9 


e OO 

5.82 


O AO 

-3.48 


A 1 A 

-0.19 


A 1 O 

0.18 


181 Ei 

hu 


1385.28 


1 OA "1 f 

1391.6 


1 O O 1 f 

1381.6 


f OO 

0.32 


O f O 

-3.68 


A 1 O 

-0.18 


A 1 C 

0.16 


182 

hu 


1 O O O 1/1 

1388.14 


1 OAO O 

1393.3 


1 O OO T 

1382.7 


E 1 

5.16 


E A A 

-5.44 


-0.16 


A 1 A 

-0.19 


183 p.. 
rjU 


1 OQ1 OO 

loyi.tjz 


loyo.o 


i qoK a 
looD.4 


^ /IS 


-o.yz 


-U.lo 


-U.lo 


171 Gd 


1362.7 


1367.3 


1365.5 


4.61 


2.81 


0.34 


0.28 


172 Gd 


1366.46 


1372.3 


1370.5 


5.85 


4.05 


0.33 


0.28 


173 Gd 


1370.19 


1375.5 


1373.1 


5.31 


2.91 


0.33 


0.28 


174 Gd 


1373.93 


1380.2 


1377.1 


6.27 


3.17 


0.32 


0.28 


175 Gd 


1377.53 


1383.1 


1379.6 


5.57 


2.07 


0.31 


0.28 


176 Gd 


1381.02 


1387.6 


1383.2 


6.58 


2.18 


0.31 


0.27 


177 Gd 


1384 


1390.3 


1385.1 


6.3 


1.1 


0.31 


0.26 


178 Gd 


1388.1 


1394.6 


1388.6 


6.5 


0.5 


0.35 


0.24 


179 Gd 


1390.11 


1397.2 


1390.2 


7.09 


0.09 


0.36 


0.22 


180 Gd 


1393.21 


1401.1 


1393.6 


7.89 


0.39 


-0.2 


0.21 


181 Gd 


1396.42 


1403.4 


1395.5 


6.98 


-0.92 


-0.2 


0.19 


182 Gd 


1399.48 


1407.1 


1398.8 


7.62 


-0.68 


-0.19 


0.16 


183 Grf 


1402.46 


1408.9 


1400.4 


6.44 


-2.06 


-0.17 


-0.19 


184 Grf 


1405.6 


1412.8 


1403.6 


7.2 


-2 


-0.15 


-0.18 


185 Gd 


1408.97 


1414.9 


1404.6 


5.93 


-4.37 


-0.13 


-0.16 



cases. This prediction not only confine to masses, radii, but also comes well for 
other observables. Thus, if we believe all these predictions as success, then the mass 
formula specially INM needed some modification, specially in the region of Z = 37 
- 40 and Z = 60 - 64 which are considered in the current work. 



6. Summary and Conclusion 

In Summary, we have calculated the binding energy, rms charge and matter radii, 
quadrupole deformation parameter for the neutron drip-line nuclei having Z = 17 
- 23, 37 - 40 and 60 - 64 regions using RMF(NL3*) which are recently predicted 
to be in islands of inversion due to their extra stability compared to the near by 
isotopes. Since the considered isotopes are experimentally unknown, we compared 
our results with various mass formula predictions. We found large differences both 
in binding energy and deformation indicating the special nature of these nuclei. 
We got some interesting features just like jerk and deep at some places in charge 
distribution radius which is different from our conventional distribution. In regions 
Z = 17 - 23, N = 42, Z = 37 - 40, N = 68, and Z= 60 - 64, N = 112 behave as more 
stable. The true properties of these nuclei can be revealed after the experimental 
observations. 
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